Abstract:
In this work, we discuss two classes of methods to reduce the complexity of (multi scale) fracture simulations. In a first part, we discuss algebraic model reduction. We show that algebraic model reduction, such as the proper orthogonal decomposition, cannot be used directly because of the lack of correlation introduced by the damage or cracks. We demonstrate the use of proper orthogonal decompositions by subdomains as a candidate to reduce computational expenses in non-linear fracture simulations whilst controlling the error level.
In a second part, we propose an adaptive hybrid multi scale method for modelling fracture in a heterogeneous material that is composed of orthotropic grains with cohesive interfaces between the grains. Instead of a direct solver, FE2 method [1] based on homogenisation is employed in order to compute the effective behaviour of the heterogeneous microscopic material on the coarser scale. At this scale, the modelling error due to the homogenisation is still low [3]. The coarse scale is discretized with unstructured triangular finite elements, and adaptive mesh refinement is used to control the discretization error. While the mesh refinement keeps the discretization error within a certain range, the modelling error increases due to the fact that by refining the coarse elements, the scale separation assumption, which is a key issue for homogenisation, may no longer be fulfilled [4].
Whereas the modelling error is inversely proportional to the size of the coarse elements, a critical element size can be found that corresponds to the critical value of the modelling error. A critical zone emerges when the size of a coarse element reaches the critical size, or if the underlying representative volume element of the microstructure loses stability due to localisation (lack of scale separation). Thereafter, a zoom-in process is triggered that replaces the corresponding coarse elements of the critical zone with high resolution micro scale mesh to which it glues the coarse scale mesh through a strong coupling technique using Lagrange multipliers [5]. The high resolution region can gradually be extended to include the newly emerging critical zones. A local arc-length technique is adopted to trace the highly non-linear curve of the global load-displacement by controlling the opening of microscopic cohesive cracks in the fully resolved regions.
The proposed adaptive multi scale method allows us to introduce progressive discrete micro cracks at the macro scale. The unstructured mesh enables us to model problems with non-regular shapes, and the arc-length method, defined over multiple scales, allows the regularisation of softening problems that are treated in quasi-statics. We exercise this method on the simulation of polycrystalline fracture, where each grain is considered orthotropic and compare results to direct numerical simulation.
|
